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Abstract 

We present an evaluation of the running coupling constant for Nf = 2 + 1 QCD. The Schrodinger 
functional scheme is used as the intermediate scheme to carry out non-perturbative running from 
the low energy region, where physical scale is introduced, to deep in the high energy perturbative 
region, where conversion to the MS scheme is safely performed. Possible systematic errors due to 
the use of perturbation theory occur only in the conversion from three-flavor to four-flavor running 
coupling constant near the charm mass threshold, where higher order terms beyond 5th order in 
the P function may not be negligible. 

For numerical simulations we adopted Iwasaki gauge action and non-perturbatively improved 
Wilson fermion action with the clover term. Seven renormalization scales are used to cover from 
low to high energy region and three lattice spacings to take the continuum limit at each scale. 

A physical scale is introduced from the previous Nf = 2 + 1 simulation of the CP-PACS/JL-QCD 
collaboration [l], which covered the up-down quark mass range heavier than ~ 500 MeV. Our 
final result is a^jg(Mz) = 0.12047(81) (48) (1?73) and A^^"^^ = 239(10) (6) (1^2 ) MeV . 
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I. INTRODUCTION 



The strong coupling constant and quark masses constitute the fundamental parameters of 
the Standard Model. It is an important task of lattice QCD to determine these parameters 
using inputs at low energy scales such as hadron masses, meson decay constants and quark 
potential quantities. The results can be compared with independent determinations from 
high energy experiments, which should provide a firm evidence of the single scale nature of 
QCD. 

In the course of evaluating these fundamental parameters we need the process of renor- 
malization in some scheme. The MS scheme is one of the most popular schemes, and hence 
one would like to evaluate the running coupling constant through input of low energy quan- 
tities on the lattice and convert it to the MS scheme. A difficulty in this process is that 
the conversion is given only in a perturbative expansion, and should be performed at high 
energy scales much larger than the QCD scale. At the same time the renormalization scale 
/i should be kept much less than the lattice spacing to reduce lattice artifacts, namely we 
require 

AQCD</i<- (I.l) 
a 

A practical difficulty of satisfying these inequalities in numerical simulations is called the 
window problem. 

One of the widely used definitions of the renormalized coupling on the lattice is to employ 
quantities related to the heavy quark potential jl, [sj , which are easy to measure accurately. 
Choosing small size Wilson loops, the running coupling constant is extracted from their 
perturbative expansion with the renormalization scale set to /i ~ 1/a. This conflicts with the 
window, and as a consequence lattice artifacts are intrinsically included in the perturbative 
expansion coefficient in terms of the MS coupling. The coefficients tend to explode, which is 
partly cured by the tadpole improvement [2!] and by combining with 0(a) improved actions 
like the staggered fermion action. This definition of the coup ling constant has been employed 
for A^f = OE, 13, = 2 0, 9 and AT/ = 2 + 1 0, i, 0, H Q flavor cases (for a review see 



Ref. [13|, uMX- Recently developed methods using moments of charm quark current-current 
correlator [l5| or vacuum polarization function [l6| are also not free from the window problem 
when applying the perturbative expansion while reducin g th e lattice artifact. 

The Schrodinger functional (SF) scheme 17, 18, fiol. 21| is designed to resolve the 



window problem. It has an advantage that systematic errors can be unambiguously con- 
trolled. A unique renormalization scale is introduced through the box size L. A wide range 
of renormalization scales can be covered by the step scaling function (SSF) technique, which 
exempts us of the requirement to satisfy the condition (II. ip in a single simulation. This 
matches our goal to obtain the coupling constant in the MS scheme and make compar- 
isons with high ener gy i nputs. The SF scheme has been applied for evaluation of the QCD 
coupling for A^/ = [Si and A^/ = 2 |2l|. 

In the SF scheme we start with the evaluation of the running coupling constant for a 
variety of the bare coupling constant (3 and box sizes, which covers the strong coupling 
region corresponding to the energy scale /x ~ 500 MeV and the weak coupling region around 
/i ~ 40 GeV. At low energy scales we expect the strange quark contribution to be important 
in addition to those of the up and down quarks. Thus the aim of the present paper is to 
go one step further than those of Refs. [H, 21] and evaluate the strong coupling constant in 



N f = 2 + 1 QCD. For setting the physical scale we employ a recent large-scale N j = 2 + 1 
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lattice QCD simulation employing non-perturbatively 0(a) improved Wilson quark action; 
the work of CP-PACS/JL-QCD Collaboration with relatively heavy pion mass with m^r ~ 
500 MeV \h]. 

II. SCHRODINGER FUNCTIONAL FORMALISM AND ACTION 



The Schrodinger functional is defined on a finite box of size x T with the Dirichlet 
boundary condition at the temporal boundary. For QCD the Dirichlet boundary condition 
is set for the spatial component of the gauge link 



f^fc(a;)Uo=o = exp {aCk] 
(<Pi 

4>3 
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(II.l) 

(11.2) 



and for the quark fields 



^(x) |^.(,=T = 0, V^(x) |^„=o = ^{x) \^g=T 



0. 



(11.3) 



Under a mild assumption it is proven that the tree level gauge effective action has a global 
minimum around a background field IT] 



V^{x) = exp {aB^{x)) , 

Bo = 0, Bk = ^ {x^C'k + (T - a;°)C, 



(11.4) 
(11.5) 



which is uniquely given by the boundary fields (]11.2|) . The fermionic mode is shown to have 
a mass gap [l9|, so that we are able to define a mass independent scheme directly in the 
chiral limit. 



In this paper we adopt the same set up (scheme) as the Alpha collaboration [18|, l2l|| for 
the boundary link (111.21) 



(11.6) 
(11.7) 
(11.8) 
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The parameter rj is used to define the renormalized coupling constant from the derivative 
of the effective action and is set to zero in the action after taking derivative with respect 
to it. The parameter v may be used to define another renormalized quantity, but we set it 
to zero when evaluating the coupling constant. We employ the choice T = L so that the 
renormalization scale is given by the box size L. 
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We adopt the renormahzation group improved gauge action of Iwasaki given by 

^. = 4 E Wo{C,gl)ReiT{l-P{C)) + ^ ^ W,{C , gl)Re ti {I - R{C)) , (II.9) 



where 5*0 and 5*1 are the sets of oriented plaquettes and rectangles. The weight factor 



is chosen to cancel the 0(a) contribution from the boundary according to [22, [23 



0/1 



cocf ((7q) Set of temporal plaquettes that just touch 
Wo(C, (y^o) = { of the boundaries, (11.10) 

Co otherwise, 

cicf{gQ) Set of temporal rectangles that have exactly 
Wi(C, (7o) = { two links on a boundary, (H-H) 

ci otherwise. 

The bulk coefficients are set to ci = —0.331, cq + 8ci = 1. The boundary improvement 
coefficients are set to the tree-level values cf = 1 and = 3/2; it is empirically known that 
they giv e better scaling behavior than the one-loop values for the Nf = [23] and Nf = 2 
case (24| . 

We used the improved Wilson fermion action with clover term 

Sf[U,^,^]=a^Y.^{Dw + mo)iJ, (11.12) 

X 

Dw = \ {l, (V^ + V;) - aV;V^) - csY,\a,,P,,. (11.13) 

The improvement coefficient csw isgiven non-perturbatively in a polynomial form for Nf = 3 
QCD with the Iwasaki action by [25| 

Csw (go) = 1 + 0.113(7^ + 0.0209(72)(?o^ + 0.0047(27)^7o', (11.14) 

which covers 1.9 < (3 < 12.0. Although 0{a) effects in the bulk is canceled by the clover 
term, there are 0{a) contributions from the boundary for the SF formalism and we need to 
add the boundary term to cancel it. 



Soia) = a' E (ct - 1) (^(^, 1)V'(5^, 1) + H^, T - l)^(x, T - 1)) . (11.15) 

X 

The coefficient is set to the one loop value given by [2^ 



Q = 1 - 0.00881(28)^^. (11.16) 
We employ the twisted periodic boundary condition in the three spatial directions, 

7/^(x + Lk) = e'^ipix), i^{x + Lk) = e-*V(x) (11.17) 
with the same 6 = 7r/5 for all spatial directions, as was used by the Alpha collaboration 
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The renormalized gauge coupling in the SF scheme is defined from the effective action 
r[V^] at the global minimum. For numerical simulation we take the derivative in terms of 
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the parameter t] introduced in the background field 0j and define the SF couphng constant 



as 
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1 dT[V, 



(11.18) 

?7=0 



g'^{L) k drj 
where 

A; = 12(^^) (co(sine + sin20+4ci(sin2e + sin40), ^ = [yl) ^"'^^^ 
is a normahzation coefficient evaluated at tree level. 



III. OUR STRATEGY 

Our goal is to derive the renormalization group invariant (RGI) scale Aqcd in physical 
units and evaluate the running coupling constant as(M^) at high energy scale /i = Mz- The 
RGI scale A is scheme dependent and we employ the commonly used definition for the SF 
scheme, 

As. ^ i m-^' exp (-^) exp (- f ^' ( . - . („,,) 

where 'g{L) is the SF renormalized coupling at the box scale L and j3{g) is the renormalization 
group (3 function in the same scheme whose perturbative expansion coefficients are given by 
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m = [bo + hg' + 62/ + ■■■], (III.2) 
^ ^11 - Inj) , (III.3) 



(47r)2 V 3 



bi = 7TIU ( - -Nj ) , (III.4) 
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,^ , , ,0.483(7) - 0.275(5)iV/ + 0.0361(5)A^; - 0.00175(1)A^? . (III.5) 
(47r)'^ ^ J J J 

The derivation of the RGI scale for the SF scheme proceeds in the following steps [jj 

(i) We start by calculating the step scaling function (SSF) T,{u,a/L) on the lattice at 
several box sizes and lattice spacings. The SSF gives the relation between the renor- 
malized coupling constants when the renormalization scale is changed by some factor, 
which is fixed to 2 in this paper, 

u,^)=t{2L) (III.6) 



L 



The scale is given by the box size L, and a/L represents the discretization error. We 
take sufficient number of values for the coupling u to cover low to high energy scales. 
Taking the continuum limit at each scale u 

aiu)= lim (III.7) 

and performing a polynomial fit we obtain a non-perturbative running of the coupling 
constant in the SF scheme for the scale change of 2. 
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(ii) In the second step we define a reference scale Lmax tlirougli a fixed value of the renor- 
malized coupling constant ^^(Lmax)- The value of ^^(Lmax) is arbitrary as long as it 
is well in low energy region to suppress lattice artifacts with a/Lmax -C 1. We then 
start from Lmax and follow the non-perturbative RG flow in the SF scheme into the 
high energy region. A typical scale turns out to be l/Lmax ~ 0.5 GeV in this paper 
so that after n ~ 5 iterations the scale 1/L = 2"/Lniax ~ 16 GeV is already in the 
perturbative region where the difference between perturbative and non-perturbative 
RG runnings is negligible. 

(iii) Substituting g'^{L) and L = 2~"'Linax into the definition (lIII.ip and evaluating the 



integral with three loops /^-function in the SF scheme [271] for the weak coupling region 



we obtain the RGI scale AsF-^^max in terms of the reference scale. 

(iv) In the last step we need some physical input measured in an independent large scale 
simulation at some lattice spacing a to quote Lmax in physical units. The requirement 
for the lattice spacing and the reference scale is that the magnitude of lattice artifacts 
O'/Lranx should bc kept small. In this paper we employ hadron masses for physical input 
and use the lattice spacing determined from them in physical units as the intermediate 
scale. We then obtain the RGI scale Asf in physical units. The transformation into 
the MS scheme is given exactly at one-loop order via 

Aj;jg = 2.61192AsF (III.8) 

for three flavors. 

The RGI scale Aj^ measured so far is for three flavors (A|^). In order to evaluate the 
coupling constant as{Mz) at high energy we need to change the number of flavors at charm 
and bottom quark mass thresholds, obtaining A^ for five flavors. For this purpose we 
used the matching formula near mass thresholds for the MS scheme at three-loop order in 
Refs. HHii. The evaluation of as{Mz) will proceed in the following steps in this paper. 



(i) Introduce the physical scale through hadron masses and evaluate Lmax in units of GeV. 

ii) Perform the non-perturbative step scaling n = 5 times and reach deep into the per- 
turbative region g ~ 16 GeV. 



(iii) Change the scheme to MS according to the two-loop relation [27 



"ms(^9) = "SF(g) + ci(s)agp(g) + C2(s)a^p(g) + ■ ■ ■ , (III.9) 
ci(s) = -Svrfoo ln(s) + 1.255621(2) + 0.0398629(2)Ar/, (III.IO) 
C2(s) = ciisf - 327r% ln(s) + 1.197(10) + 0.140(6)A^/ - 0.0330(2)Ar/. 

(III.ll) 

We may set the scale boost factor s = 2.61192 so that ci{s) = 0. A systematic error 
due to higher loops correction is less than 0.1 % and negligible here. 

(iv) Running back to the charm quark mass threshold /i = rric with the four loop jS- 
function in the MS scheme we change the number of flavors to four using the three-loop 
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matching formula 28|, |29|, |30 



TT TT ^ ^ jl^ 

3 /Q;\ 

= l + Y.Fu{x)(-] , (III.13) 

fc=i 

i^i(x) = ^x, (III.14) 
o 

F^ix) = Fi(a;)2 + llx + ^, (III.15) 

564731 82043^,^, 955 53 2 1 ^ 

FJx) = C(3) H X H H a;^ 

^ ^ 124416 27648^ ^ 576 576 216 

where M{ii) is the MS running mass of the heavy quark which decouples at the thresh- 
old. We shall set /i = M(M) and x = in this paper. Since the largest error may 
be introduced from the use of perturbation theory at /i = mc(mc), we estimate the 
systematic error of this perturbative matching, by comparing the result with that from 
the two-loop matching relation [sil, 32. 33l . 

(v) Running to the bottom quark mass threshold /x = mf,(mf,) we obtain the running 
coupling constant for five flavors in the same manner. 

(vi) Finally we change the scale to /i = Mz{Mz) with the four-loop /5-function and find 
Ois{Mz). 

vii) The RGI scale is given by substituting /i = Mz{Mz) = 1/L and as{Mz) in the 
definition fllll.ip for five flavors in the MS scheme with the four-loop (3{g). 



IV. STEP SCALING FUNCTION 



We adopt seven renormalized coupling values to cover weak = 1.001) to strong = 
3.418) coupling regions, which approximately satisfy gfj^^{L) = g1{2L) (i = l,---,6). For 
each coupling we use three boxes L/a = 4, 6, 8 to take the continuum limit. 

The HMC algorithm is adopted for two flavors and the RHMC algorithm for the third 
flavor, all of which are set to a common mass of zero. We adopt the CPS++ code [34] and 
add some modification for the SF formalism. Simulations were carried out on a number of 
computers, the PC cluster Kaede, PACS-CS and T2K-tsukuba at University of Tsukuba, 
T2K-tokyo and SRllOOO at University of Tokyo and the PC cluster RSCC at RIKEN. 

The distribution of the inverse of the coupling constant 1/^^ turned out to be a smooth 
Gaussian even at the lowest energy scale |2J] as plotted in Fig. [H This is contrary to the 



finding with the standard Wilson gauge action [18|, |21j and we need no re-weighting. 

We start by tuning the value of f3 and k to reproduce the same renormalized coupling 
at each of the box sizes 4, 6, 8 keeping the PCAC mass to zero. Requirement for the 
renormalized couplings g^{L) is that their values agree within one standard deviation for 
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L/a = 4, 6, 8. The PCAC relation is defined in terms of the improved axial current with 
non-perturbative improvement coefficient 351 



A'rix) = A,{x) + CAd,P{x), CA{gl) = -0.0038^g |_°-^^^J . (IV.l) 

The values of k) are listed in Table [T] together with results for the renormalized coupling 
constant ^ and the PCAC mass at the two scales L and 2L. Statistics of the runs are given 
in Table M 

The renormalized coupling ^^(2L) at the scale 2L is corrected perturbatively in order to 
cancel the deviation of the PCAC mass from zero at the scale L[36| 

^'(2^)-2..._ _= ^'(2^) -$(OymL, (IV.2) 



g (L)=u,m=0 



g {L)=u,m{L)=m 



$(0) = 0.00957A^/. (IV.3) 

The PCAC mass at the scale L has been tuned such that the deviation ^{0)u^mL is smaller 
than the typical statistical error. 

The value of the renormalized coupling g'^{L) at L/a = 8 is used to define g'^{L) at scale 
L. The deviation of g'^{L) at L/a = 4, 6 from it is also corrected perturbatively at three- loop 
using 
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Lj \ Lj du 

so = 26oln2, (IV.6) 
si = (26o In 2f + 2hi In 2, (IV. 7) 

S2 = (26o In 2)^ + lOfoo&i (In 2)^ + 2^2 In 2, (IV.8) 

where bn is the perturbative coefficient of the /3-function in the SF scheme. 

We now consider the continuum extrapolation a/L of the SSF. In perturbation 
theory the deviation of the lattice SSF from its continuum value is expressed as 

\\ — = 5i{a L)u + 52ia L)u^ + ■ ■ ■ , IV.9 

a{u) 

5i{a/L) = 5iG{a/L) + Nf5iQ{a/L). (IV.IO) 

The one-loop coefficients 5ig/iq are given in Table UTTl for the Iwasaki gauge action with the 
tree-level improved boundary coefficients q adopted for the present work for each box sizes 
2^, 13] • As is seen from the table the values of (5iq/ig are not small, and the deviation 



decreases only slowly with the volume L/a. 

Instead of calculating the two- loop coefficients ^2QI2G perturbatively, which is a non- 
negligible task, we calculate SSF directly by Monte-Carlo sampling at very weak coupling 
/? > 10. The results are listed in Table llVt where the parameter is tuned only for k to 
reproduce mpcAc = 0. We define the deviation from the perturbative SSF 



ll{u,a/ L) — crp^(-u) 
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where ap^ is the continuum SSF at three- loop order given by fllV.Sp . The deviation is fitted 
in a polynomial form for each a/L, 

1 + 5{u, a/L) = 1 + di{a/L)u + d2{a/L)u^. (IV.12) 

We tried a quadratic fit using data at m < 1.524 with fixing di{a/L) to its perturbative 
value 6i{a/L), which is plotted in Fig. [2l We also plot perturbative one loop behavior for 
comparison. As is seen from the figure the one loop line could reproduce the data only at 
very high /3 > 10 for L/a = 4, 6. It may not be safe to adopt the one loop improvement for 
our data at m > 1.0. 

The fit results for the coefficients are listed in table El We observe that the higher-loop 
coefficient ^2 is not negligible and contribute in opposite sign. We notice that the fit result 
hardly changes even if we add one more data at u = 1.840. Since the quadratic fit provides 
a reasonable description of data as shown in Fig [2] we opt to cancel the 0{a) contribution 
dividing out the SSF by the quadratic fit according to 

E(2) (u -) = ^(^^(^/L) . 

V'LJ l + 6i{a/L)u + d2{a/L)u^' ^ 

Now we have the values of the 0{a) improved SSF's in the chiral limit for three lattice 
spacings at each of the 7 renormalization scale given by u, which are listed in table IVIi 
Scaling behavior of the SSF is plotted in Fig. [3l Almost no scaling violation is found. We 
performed three types of continuum extrapolation: a constant extrapolation with the finest 
two (filled symbols) or all three data points (open symbols), or a linear extrapolation with 
all three data points (open circles). As is shown in the figure they are consistent with each 
other. Since the scaling behavior is very good for the finest two lattice spacings we employed 
the constant fit with these two data point to find our continuum value, which is also listed 
in Table IVD 

The RG running of the continuum SSF is plotted in Fig. HI We divide the SSF with 
the coupling g'^{L) to obtain a better resolution in this figure. A polynomial fit of the 
continuum SSF to sixth order fixing the first and second coefficients Sq and Si to their 
perturbative values flIV.6p . flIV.7p yields 

a{u) = U + Squ"^ + SiU^ + S2U^ + SsM^ + 54^*^, (IV.14) 

52 = 0.002265, S3 = -0.00158, S4 = 0.000516. (IV.15) 

The fitting function is also plotted (solid line) together with the three loop perturbative 
running (dashed line). 



A. Non-perturbative /3-function 

From the polynomial form of the SSF we derive the non-perturbative /3-function for 
Nf = 3 QCD. Starting from definition of the /3-function 

-L^^ = 2V^(3{V^), u = f{L) (IV.16) 

the value of the /3-function at stronger coupling (lower scale) is given by recursively solving 
the relation 
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The input is the three loops perturbative value at u = 0.9381, which is deep in the pertur- 
bative region. 

For the non-perturbative SSF we adopt a slightly different fitting form in order to reduce 
the error propagation. We performed a polynomial fit by fixing the first to third coefficients 
So, Si and S2 to their perturbative values (lIV.6p . (lIV.7p . (IIV.SP 

a{u) = U + SqU^ + SiU^ + S2U'^ + S^U^ + S4U^, (IV. 18) 

S3 = -0.000673, S4 = 0.0003434. (IV.19) 

The resultant /^-function is plotted in Fig. [51 The /3-function oiNj = 2 QCD is reproduced 
from data of the Alpha collaboration [2]]] for comparison. Note that the error is estimated 
by a propagation from those in the continuum SSF's a{u). 



V. INTRODUCTION OF PHYSICAL SCALE 

CP-PACS and JLQCD Collaborations jointly performed an Nf = 2 + 1 simulation with 
the 0(a) improved Wilson action and the Iwasaki gauge action, whose results have been 
recently published [l[. Three values of /5, 1.83, 1.90 and 2.05 were adopted to take the 
continuum limit and the up-down quark mass covered a rather heavy region corresponding 
to m^/mp = 0.63 — 0.78. 

We adopt those results to introduce the physical scale into the present work so that the 



reference scale Lmax is translated into MeV units. The Alpha Collaboration [18|, |21| has 
adopted the Sommer scale tq as a physical observable for this purpose. Since the Sommer 
scale is not a direct hadronic observable, we prefer to employ the hadron masses m^, m^, 
as inputs and use the lattice spacing a as an intermediate scale, which are listed in Table 

We evaluate the renormalized coupling in the SF scheme at the same (3 = 1.83, 1.90, 
2.05 in the chiral limit. The reference scale Lmax is given by the box size we adopt in this 
evaluation. Note that this definition gives a different value of //max at different (3. The 
renormalized coupling ^^(I/max) should not exceed our maximal value 5.13 of the SSF very 
much. The value of the coupling constant at each f3 are listed in Table IVIIII together with 
the PCAC mass. The hopping parameter n is tuned to reproduce mpcAc = except for the 
cases that the couphng constant apparently exceeds 5.13. We use the box size of L/a = 4 
for (3 = 1.83 and 1.90 to define L^ax and L/a = 4, 6 for /5 = 2.05. 



VI. RGI SCALE AND THE STRONG COUPLING CONSTANT AT Mz 

Starting from Mmax = ^^(/^max) we iterate the non-perturbative renormalization group flow 
five times according to the polynomial fit flIV.14|) and substitute the result L = 2~^Linax 
and g{L) into fllll.ip with /3-function for three flavors at three loops. In this way we obtain 
AgpLinax for three flavors. Further non-perturbative step scaling with n > Q does not change 
the central value of AgpLmax- The results are listed in Table HXl together with Agp in units 
of MeV and given by fUTLSD . 

We derive the strong coupling constant as\Mz) at high energy scale ^ = Mz according 
to the procedure given in Sec. Illli After reaching the scale L = 2~^Linax in the SF scheme, 
we transform to the MS scheme by the two-loop formula (1III.9P at q = 1/L with s = 
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exp(ci(l)/(87r6o))- Then running back to the scale /i = mc{mc) with three-flavor 4-loop (3- 
function the coupling constant is matched to that for four flavors at three-loop order using 
(]III.12p . We repeat the same operation at the threshold n = mb{mb) and obtain the five flavor 
coupling constant. We finally run to fi = Mz with the four-loop /3- function for five flavors 
and find as(M^). The QCD parameter is given by substituting /i = Mz = 1/L and 
as{Mz) in ( jlll.ll) for the MS scheme with 4-loop l3{g). The results are listed in Table 1X1 For 
an estimate of the systematic error due to perturbation theory, results using three- and two- 
loop formula in (]III.12p are listed. The error includes the statistical error of the renormalized 
couplings, which is propagated into that of the SSF, in addition to the statistical error of 
the lattice spacing. The experimental errors of rric, mi, and Mz are also included. 

As the last step we take the continuum limit using the three lattice spacings from Ref. 
The scaling behavior of as{Mz) and is plotted in Fig. [61 Since the results in the 
continuum limit do not depend on Lmax, we adopt the result for L = 6 as the central value 
for /5 = 2.05. 

We tested three types of continuum extrapolation; a constant fit with three or two data 
points, or a linear extrapolation ^. These results agree with each other and we adopt the 
constant fit with three data points for our final results since there is almost no scaling 
violation. Our final results are 



where the first parenthesis is statistical error and the second is systematic error of per- 
turbative matching of different flavors, which is estimated as a difference between results 
with three- and two- loop matching relation for (1111.120 and may be overestimated. The 
last parenthesis is a difference between the constant and a linear extrapolation and is a 
systematic error due to finite lattice spacing for physical inputs. 

VII. CONCLUSION 

We have presented a calculation of the running coupling constant for the Nf = 2 + 1 QCD 
in the mass independent Schrodinger functional scheme in the chiral limit. We used seven 
scales to cover low to high energy regions and three lattice spacings to take the continuum 
limit at each scale. 

After tuning j3 and /t to fix seven scales in the massless limit we evaluated the step scaling 
function in the continuum limit. We notice that deviation fllV.lOl) from the continuum SSF 
is rather large at one loop for our choice of the Iwasaki gauge action and the tree level 
improvement for boundary coefficient c^^^. Since the one loop formula could not reproduce 
the numerical data except for very high /3 > 10 we adopted "two loops" formula extracted 
from numerical data with quadratic fit. With the "perturbative" improvement the SSF shows 
good scaling behavior and the continuum limit seems to be taken safely with a constant 
extrapolation of the finest two lattice spacings. 



0{gQa/L) error is expected from boundary terms in temporal direction in the SF scheme, which may 
propagate to as{Mz) through 5^(Linax)- 




(VI.l) 
(VI.2) 




239(10)(6)(1°2) MeV, 
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We notice that "two loop" term in the deviation (1IV.11|) has been comparable to that at 
one loop. There may be a possibility that higher order perturbative correction contribute in 
an non-negligible manner, which may introduce an unestimated systematic error. However 
we consider the probability is not so high since scaling behavior of the "two loops" improved 
SSF is good as in FigJH] and the continuum limit was taken safely. But a further test may 
be preferable with a different setup with better perturbative behavior for the SSF. 

With the non-perturbative renormalization group flow we are able to estimate the renor- 
malization group invariant scale Aqcd and as{Mz) with some physical inputs for energy 
scale. The physical scale is introduced from the spectrum simulations of CP-PACS/JLQCD 
collaboration [T] through the hadron masses m^r, mx, ttlq. From these inputs we evaluated 
(IVI.ll) and flVI.2p . where all the statistical and systematic errors are included. Our result is 
consistent with recent lattice results [l^, 11, 12, 15 1 and the Particle Data Group average 



as{Mz) = 0.1176(20) [38|] with the systematic error included. 

For a future plan a new result is going to be available by the PACS-CS Collaboration 
i, 13] aiming at simulations at the physical light quark masses down to mT^/nip ^ 0.2. 



This may reveal a systematic error from the physical scale input due to chiral extrapolation 
toward light quark masses. With progress in the physical point simulation expected in the 
near future, we are hopeful that a full control of errors in the lattice QCD determination of 
the strong coupling constant is in sight. 
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TABLE I: The value of /? and k to reproduce the same physical box size L and near zero PCAC 
mass. The renormalized coupling and PCAC mass at scale L and 2L is also listed. 
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TABLE II: Number of configurations for each run. 
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p IP 

TABLE III: Perturbative improvement factor at one loop level for tree level improved Cj . Pure 
gauge contribution b\G is taken from Ref. 20] and quark contribution b\Q from 37 1. 
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TABLE IV: The renormalized coupling and PCAC mass at scale L and 2L to derive the SSF at 
high P > 10. 
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TABLE V: Coefficients of the quadratic fit of the deviation 6{u,a/L). 
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TABLE VL The 0{a) improved SSF 12^'^\u,a/L) at "two loop" level for three lattice spacings 
a/L = 1/4, 1/6, 1/8. Corrections are made such that the PACS mass m = and the renormalized 
coupling u at smaller box to be the same value for each lattice spacings. The SSF a{u) in the 
continuum is also listed, which is given by a constant fit of two data at finest lattice spacings 1/6, 
1/8. 



17 





1.83 


1.90 


2.05 


a (fm) 


0.1209(16) 


0.0982(19) 


0.0685(26) 



TABLE VII: Lattice spacing a from large scale simulation [l| 
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TABLE VIII: The renormalized coupling and PCAC mass at /3 = 1.83, 1.90, 2.05 used in the large 
scale simulation. The values of k's are tuned to reproduce mpcAC = except for the case that the 
coupling apparently exceeds 5.13. 
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TABLE IX: The RGI scale for three flavors in the SF scheme and Ai^ in the MS scheme. 
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TABLE X: The strong coupling as{Mz) at /i = Mz and the RGI scale Ai^^ for five flavors. Those 
in the third and the fourth column are derived with three loops formula for piI.12p . Those in the 
fifth and the sixth are from two loops formula. 
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FIG. 1: Distribution of inverse of the renormalized coupling at lowest energy scale given by g'^{L) ~ 
5, which corresponds to L/a = 16, /? = 2.5 (left), L/a = 12, /? = 2.34652 (middle) and L/a = 8, 
P = 2.15743 (right). Solid line is a fit in a Gaussian function. 




PT one loop 
- quadratic fit with fixed d, 



PT one loop 
- quadratic fit witfi fixed d, 



0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 



FIG. 2: Polynomial fit of discrepancy S {u,a/L) /a"pr|(rt) at high (3^4. The fit is given for each 
lattice spacings a/L = 1/4 (left), a/L = 1/6 (middle) and a/L = 1/8 (right). Black dotted line is 
a perturbative one loop behavior and red solid line is a quadratic fit with fixed di to its one loop 
value. 
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FIG. 3: The SSF on the lattice with its continuum extrapolation at each renormalization scale. 
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FIG. 4: RG flow of the SSF divided by the coupling g'^{L). Dotted line is three loops perturbative 
running. Solid line is a polynomial fit of the SSF. 
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FIG. 5: Non-perturbative /3-function for Nf = 3 and 2 QCD. Solid lines are three loops perturbative 
running for comparison. Data for Nf = 2 is reproduced from Ref. 21 1. 
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FIG. 6: Scaling behavior of ay^{Mz) (left) and (right). We adopt 6^ data for /3 = 2.05. 
Three types of continuum extrapolation is given; constant fit with three and two lattice spacings 
and linear fit. 
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